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CN . Abstract 

t~^ ■ The connection between modulated Riemann surface of genus one 

^N , and solution to Volterra lattice that tends to constants at infinity is 

studied. The main term of asymptotics for large time of solution to 

the mentioned Cauchy problem is written out. 

c/5 '. 1 Introduction 
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The Volterra lattice 

DtCn = Cn{Cn+l - C„_i), C„ = C„(t), 71 G Z (1.1) 

is known as an interesting example of integrable difference- differen- 
tial model with various applications to plasma physics, crystallogra- 
phy, zoology. The Inverse Scattering Data Method (ISDM) allows to 
investigate in details Cauchy problem for the Volterra lattice in case 
of quickly decreasing initial data (see ^) and in periodical case (IQ], 
[^ ) . The further interest is stimulated by problems with more complex 
behavior of solutions to (|1.1[) as n ^ ±00. 



Examine the Cauchy problem for eq (LI) with initial data (i.d 



Cn(0) which quickly tend to constant values as n ^ ±00: 



"?<°'='-<°> = "',. "-±» (1.2) 

[ Vn{0) = C2„+l(0) =V^ j 

The goal of this paper is to construct asymptotics as t ^ oo for 
solution of problem ( |l.l[ )-(pr^). Some similar problems were studied 
for Korteweg-de Vries equation (KdV) (0]), Nonlinear Schrodinger 



equation ([|]), Modified KdV (@). The problem (|0|)-(|Tj) was posed 
in where the author formulated hypothesis to be proved lower. 

The object of study of the next paragraph is a linear equation 
associated with Volterra lattice together with appropriate initial data 
for using the ISDM. In the third paragraph we formulate Riemann 
problem corresponding to initial conditions ( |1.2| ) as u^ = v^ = c^. 
The last restriction makes calculations much less tedious whereas all 
the results can be easily expanded to the general case u^ ^ v^. 

In the fourth paragraph we describe one-gap solutions to the Volterra 
lattice, modulated Riemann surface and its connection with the one- 
gap solutions. 

The last paragraph contains asymptotic as i ^ oo investigation 
of the formulated above Riemann problem. We prove a statement on 
representation of the main term of asymptotic solution to the Cauchy 
problem ( |1.1[ )-(|L^) by modulated one- gap solution to the Volterra 
lattice. 

2 Inverse Scattering Data Method 

Examine a linear problem 

cy_\V'n-l+cy^V'n+l = AV'n, Cn -^ C±, n ^ iboo, C+ > C- > 0, (2.1) 

Cn are solutions to Volterra lattice, A is a spectral parameter. Define 
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Riemann surface of genus r-|-(2:-|-. A) : 2Ac_|_ = z±+z^ . Involution 
cr{z±) = z^ near circle |^;| = 1 induces change of sheet of the surface 
T± near the cut [0,c-i-] = E±. 



Fix solutions of problem (2T) by asymptotics 



^tW ~ et{z±) = z^'^exp Qte± (4 - z±2)^ , n ^ ±oo, (2.2) 



tp^{P) is analytic on the lower sheet T^ \oo~ while ip~{P) is analytic 
on the upper sheet F^ \ oo"*". In the spectral domain £'_)_ there is a 
scattering correlation 

V'+(P) = a(P)V-(P) +6(P)^-(^P), P e E+. (2.3) 

Quantities a{P), b{P) are connected by condition 



(2.4) 
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W {'ijjnift^n) = Cn (V'n0n+i-V'n+i0n) IS & Wronskian. There is also 
equality 



"v^'' - 1727 ITT ' ^ ' ~ Wt ^T^' ^ ' 

c_ (z_ — 2;_ j c_ (z- — z_ j 

which signifies that function a{P) can be analytically expanded to the 
lower sheet r~. 

We suppose the spectrum to be solitonless, i.e. a{P) has no zeroes 
on r~. One easily verifies that under P G E^, = E^\E- the condition 
a{P) = b{aP) holds whence (as well as from (p.4[)) one obtains: 



r{P)r{aP) = I, P £ E^, (2.6) 

where r{P) = h{P)/a{P) is reflection coefficient. 

3 Riemann problem. 



There is a general scheme for investigation of integrable systems via 
solving the inverse scattering problem in the form of matrix Riemann 
problem (MRP). 

1. Let ^n{P) be 2*2-matrix piecewise analytic in C function with 
asymptotics under P -^ oo-t 



%,(P) =A±{I + X-'V^ + 0{X 



where 



1-2^ 



(Ac±^/^)^" 



Ac 



-1/2 



ibn. I ) 



(3.1) 



A -^ oo-t 



1 



A± 



2. Let function \I'„(P) be analytic outside of the spectrum E-^- and 
on cut E^ have a jump 

*„(P - iO) = ^n{P + iO)G{P), P G E+, where (3.2) 



G(P) 



r-(P) 

r(P)r(o-P)-l 

A- 



r-(c7P) 
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, A = c_'"(z--zZ^), (3.3) 



the upper index (-) signifies value on lower side of cut £'+. It fohows 
from ( ^.6| ) that matrix G{P) becomes triangular in spectral branch 

3. Non-degenerate matrix ^„(P) is a regular function in C. 

Statement 3.1. Solution c„ of Cauchy problem ( p^ ) - (L2) for 
Volterra lattice is determined by solution ^„(P) of formulated above 
MRP (|1D- (U) via the following formula: 



K 



n-l 



21 



Vr^ 



21 



(3.4) 



where T/"*" is matrix coefficient of expansion ( 3.1 



Proof. One must just verify that exact solution to MRP ( |3.1| ) - 
( |3.2D is given by matrix 



^n(P) 



V'+(P) V'^-(aP) 

^n+l(^) V'-+l(^^) 



vF-i(^^(^),V';r(^^)) o\ 

l) 

(3.5) 



Statement 3.2. Solution of the MRP ( p.l| ) - ( |3.2D is unique. 

Proof. Let * and ^ be two solutions of MRP (^!ll) - (|]|). Accord- 
ing to Liouville theorem the matrix ^(P)^~^(P) does not depend on 
P. Thus, the quantity V^^ is derived from the MRP solution uniquely. 

Investigation of solution to Volterra lattice as t ^ oo implies 
asymptotic analysis of MRP (O) - (PI). 



4 One-gap solutions of Volterra lat- 
tice and modulation equations 

Exact formulas for real solutions to Volterra lattice were written out 
in paper |^]: 

C2n = Un = u{t) = C(2tjr) - C{2uT - a_) - C{a+) - ({a- - a+), 
C2n+i =Vn = v{t) = C{2ujT - 0+) - ClSu;'?") + C(a-) - C(a- - a+), 

(4.1) 
where r = ^{a- — a^)n -\-t + uj'; ({x) is Weierstrass zeta-function: 
C'{x) = -p{x), X = /^(^) -y=^ ; p{x + 2co) = p(x); p{x + 

2a;') = p{x); a±, lo, ioj' are four real parameters completely deter- 



mining formulas (4.1). Define the following four quantities: 



r, =2C Q(a_ -a+) +cUj)-2?7,+C(a+)-C(a_), j = 1,2,3,4, (4.2) 

where wi = 0, u;2 = w, uj-^ = uj + uj' , uja = uj' , rjj = Ci^j)- 

The four parameters r = (ri,r2,r3,r4) also completely specifies 

solutions (41) u„ = M(r, r), f „ = v{T,r). 

Suppose quantities r depend on variables x S R, t : r = r(x, t). 
Definition. On Riemann surface T{w,X) : nP' = Ri{X) = (A — 

ri)(A — r2){\ — r^){\ — r4) we define a couple of Abel differentials 



^Q='^{\ + dQ)w-^d\ ni = ]^(\^-]^k\ + d^w-^d\, (4.3) 
where A = ri + r2 + ra + r4; do, di are constants fixed by conditions 



/ no= Jii = 

Jr2 Jr2 



(4.4) 



Equation 



DtQn = DM 



x^n 



(4.5) 



is called Whitham equation or modulation equation. 

Formula (|4.5| ) can be rewritten in the form of system of quasilinear 
differential equations for which the quantities r-,-, j = 1,2,3,4 are 
Riemann invariants: 



Dtrj = Wjir)D,rj, j = 1, 2, 3, 4, (4.6) 

where Wj{r) = Qi/QQ{rj). In paper Q characteristic velocities Wj 
were written out explicitly in terms of complete elliptic integrals. In 
the same paper self-similar solutions to system ( |4.6| ) were computed. 
Here we are interested in solutions that depend only on unique variable 
^ = x/t : r = r(^). 

Define Riemann surface T(w, A; ^) in the following way: 



Tiw,X;0 



w' = iX- ri)(A - r^){X - r^{0)i^ - u), T < ^ < Co 

w^ = ix- n){x - r^{o)H^ - n), eo~ < e < eo+, 

(A - ri)(A - 4{0){X - 4){X - u), ^0+ < e < e+ 



w 



w' = {X-r^){X-u), ^>e, 



(4.7) 



where ri < 7^2 (0 < r^ < ^2(0 — '''2 — ^3 iO — '^4! ^1)^3 ) ''^2 ''^4 



are 



constants, 



{u-ri){r4 -r^) 
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log 



yjrj-ri + yjrj-r^ 



v'r4-ri--\/r4-r 



(r^ -ri)(r4 -ri) r^ - n 



(^+ + r4 - 2ri) - 2 ^ - — ' ^ + 
I ^4 + To — 2ri 



log 



^0 = 2^^i + ^-^^ ~ ^^2 > ^o' = ^{n + r^) - r^ 



]_r4-r£ 
2 7-4 -ri 



Dependence of quantities r^ (0 ) ^^2 (0 1 '''2' (0 ^^ variable ^ is deter- 
mined by self-similar solutions to system ( ^4.51) : 



e + W^3(ri,r2-,r3-(0,r4) = 0, T < ? < ^o", 
^ + W3{n,4iO,r^2iO,rA) = 0, eo" < e < ^o"", (4-8) 

^ + W2iri,4{0,rt,r4) = 0, ^0+ < C < ^+ ■ 
Quantities r = (ri, r2, rs, r4) as functions of variable ^ are depicted 
in Fig.l. 

The Riemann curve combined in such a way is a generalization of 
that obtained in paper |0]. Define basic cycles as it is shown in Fig. 2. 
Holomorphic differential Q = Dw^^dX ,D is a constant fixed by 
condition <f^Q = 1. 

Examine the following ansatz for solutions of Volterra lattice: 

c{t I r{w, A; 0)x=n + Oir^), 6 > 0, (4.9) 

( u(T,r(^)) 

f(^) are branch points of curve V{w^\\^. Quantities u^ ^v^ ^r are 
connected in the following way: 





4ri = u + w — 2\/utj, Avo = u + v + 2v' 



u V 



(4.10) 



4r^ = ■U+ + t)+ - 2\/n+i;+, 4r4 = n+ -h f + + 2Vu+v+, 

whereas one suppose that ri < r^ ■ 

Investigations in this and the next paragraphs deal with case u^ = 
v^ = c^, ri = rg" = 0, a+ = — a_ which does not change principal 
scheme but considerably simplifies calculations. 

5 Asymptotic solving the Riemann prob- 
lem. 



In paper Q explicit formulas for the case of one-gap spectrum of linear 
problem (2.1)were obtained. These formulas represent Baker- Akhiezer 
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function in terms of elliptic functions: 



7^1 



^{t, z) = 7±^ exp [±Tri {tV{z) + nP{z))] 



9{oj't ±z + P) 



9{±z + P)9{uj'T-a + P)' 
^,(5.1) 

where a = a+ = — a_ , P = uj+a, 9{t) = (t{t+u)) exp f —tt] — ^^^ ))''"■ 

an + Vt, ^& = C (''")) variable z is connected with parameter A via 
uniformization 



X{z) = Ciz - a) - Ciz + a) + 2C(a) 



rXiz) nX{z) 

V{z) = / Qi, V{z) = / Qa 

Jr\ Jri 



(5.2) 



7n=/3" 



a{u;'T + uj)a{2uj'a + w) 
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Examine function 

■ a-i(P)v+(P), Per 



/3 = o-(2a;'a + u;)/l^ 



^n(P) 



V'-(P), PGr+. 



(5.3) 



It is analytic on each sheet of the curve and has a jump on contour 
c?r+ = £'+ equal to r{P)ijj:^{aP). So the function ^n{P) can be 
restored by formula 

^n{P) = e~{P) + ^ f M{P,Q)r{Q)^-{aQ), (5.4) 

M(P, Q) is Cauchy kernel: 

W{e-{P),e-{aQ)) d\{Q) 



M{P,Q) 



WieniQ),en{aQ)) X{P) - X{Q)' 



(5.5) 



Suppose now that ^ G 



^ , i^Q and dynamics of the branch point 



Tg (^) is unknown. Then function <1>„(P) is specified by explicit for- 
mula via Cauchy integral on curve T 



c^„(P) = e-{P) + ^ /- M{P, Q)f{Q), (5.6) 



where dT~^ 



ri,r2 



U 



^3 iO,u] ; L = [r2 , rg (O] Ucr {[r^ (C), r2 ]) , a 



is an involution that changes sheet on F, 

W{e-{P),e-{aQ)) dX{Q) 



M{P,Q) 



W{e^{Q),e;^{aQ)) \[P) - X{Q) 



(5.7) 



Choice of function / is inspired by analogy with Zakharov-Manakov 
ansatz from paper |^]: 

^-{P) = A(P)e-{P) + B{P)e-[aP). (5.8) 

Substitute this ansatz into the linear equation and get: 



A{P)e-{P) + B{P)e-{aP) = 

(5.9) 
Examine equation (|5.9|) as n — > oo and use the known formula (t — > oo) 



1 



/(^) 



27ri J-oo 1] — fj, + iO 



exp(it7(r/))d7? = -f{fi)e-Kp{itj{fi))H{--f'{fi))+0 



ViifJ- - fJ'Cr) I ' 



where 7'(/icT) = 0, i:/^(x) = (1 + signx)/2 is Heaviside function. The 
role of 7(//) is played by differential Q.q{P)^ — Q.\{P). 
Therefore one obtains: 



B{P) = -r{P)A{aP)H{\{P) - \{Pct)), P G dT+ , (5.10) 

where point Pqt is determined by condition 

^o{PcT)i = ^i{Pct). (5.11) 

Examine now ( |5.9D as n — > — oo and get integral equations on A{P): 

A{P)e-{P) = e-{P) + ^J~ M{P^Q)r{Q)r{aQ)A{Q)e-{Q)H{\{P) - X{Pct)), 



p e 5r+, 



(5.12) 



1 



A(P)e;(P) = e-(P)+— I M{P,aQ)r{aQ)A{Q)e-{Q),PeL 

2.m Jdr+uL 

(5.13) 



Thus the function / in ([5.6D must look in the following way: 



-r{P) [A{aP)e~iaP) + BiaP)e~{P)] , P e dT^ 



f{P) = \ {l + r+{aP))A{P)e-{P), Pe 

il-r+{P))A{P)e-iP), PGa{ 



^2 , ^3 (0 
r2,r3iO 



(5.14) 

The quantity A{P) should be searched from eqs (|5.12| ) - ( ^.131) . They 
indicate that function A(P) is analytic on sheets T± and has a jump 
on contour dT^ U L: 



A~{P) = AUP)9{P), PGdT-^UL, 
where conjugation matrix g{P) has the form: 



(5.15) 



1 - r{P)r{aP))H{XiP) - X{Pct)), P G df+, 
g{P) = \ -r+{aP), P G [r^,r^{0] , 

r+iP), PGcT{[r,,r,{0]). 

Solution to such a problem is given by the following formula: 



(5.16) 



A{P) = R{P) expi-— I M{P, Q) log g{Q) 

V 27ri Jar+uL 

where M{P, Q) is an analogue of Cauchy kernel: 



MiP,Q) 



miP',Q), 



OO-I- 



m{P, Q) is a meromorphic bidifferential fixed by conditions 



(5.17) 



(5.18) 



i m{P, Q) = 0, m{P, Q) 

J a 



dq{P)dq{Q) 
{q{P)-q{Q)y- 



■,P^Q, (5.19) 



g is a local parameter. It ensues from norm conditions ( 5.1S| ) that 
§f,m{P,Q) = 27rin{P). 

Multiplier R{P) is specified by condition of uniqueness of A{P) 

1/2 

and asymptotics A{P) — > c_ , P — > oo-i- : 
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R(P) = const ; „+ <-, (5.20) 

{SL, ^ + Do' 



where Dq = — ( J^ + J^ 1 Q—K,KisaL vector of Riemann constants 
(see §). 



The conjugation conditions ( 5.16 ) provide a formula for phase shift: 



D = Da-^.L Vt{Q)\ogg{Q). (5.21) 

So we constructed function ^n(-P) (|3.5D which in case r^{^) = 
X{Pct) asymptotically satisfy the conjugation condition ( |3.2D and es- 
timate (|3.l| ). 

Statement 5.1. Under r^ (^) = X(Pct) the following condition holds: 

^„(P + iO)G(P)^„(P-iO) = / + 5(P,C,t), 
where small parameter 5 has the form 

up^.^jO {t-'{P - Pct)) , \P - Pct\ > t-S 



Analogously one studies the Riemann problem ( |3.2| ) for all values of 
parameter ^. 

Now use formula (|3.4D and obtain 

Statement 5.2. The main term of asymptotics as t — > oo of solution 



to Cauchy problem (LI) - (|1.2D for Volterra lattice under condition 
of solitonless spectrum is described by Whitham-modulated one-gap 
solution c{n,t \ r{^),D{^)) ([4.lD , (|4.2|) , where Riemann curve r(,^) is 
determined by formulas ( [4.7D , phase shift -D(^) - by formula ( |5.21| ). 
The remaining terms of the asymptotic series decrease as powers in t. 
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